Linear entropy in quantum phase space 
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We calculate the quantum Renyi entropy in a phase space representation for either fermions or bosons. This 
can also be used to calculate purity and fidelity, or the entanglement between two systems. We show that it is 
possible to calculate the entropy from sampled phase space distributions in normally ordered representations, 
although this is not possible for all quantum states. We give an example of the use of this method in an exactly 
soluble thermal case. The quantum entropy cannot be calculated at all using sampling methods in classical 
symmetric (Wigner) or antinormally ordered (Husimi) phase spaces, due to inner product divergences. The 
preferred method is to use generalized Gaussian phase space methods, which utilize a distribution over stochastic 
Green's functions. We illustrate this approach by calculating the reduced entropy and entanglement of bosonic 
or fermionic modes coupled to a time-evolving, non-Markovian reservoir. 



I. INTRODUCTION 

Quantum dynamics and thermal equilibrium states in large 
many-body systems have been widely investigated using 
phase space representationsll|-|4(]. Applications of these 
methods[5] include EPR correlations in parametric ampli- 
fiers, quantum soliton propagation[6], non-equilibrium quan- 
tum criticality H, quantum dynamics of Bose-Einstein con- 
densates in two[8] and three dimensions!^, flOll . molecular 
downconversion ll ill and many other problems. Gaussian op- 
erator expansions llL2T - [l5ll are an important extension to these 
phase space mappings. Unlike earlier methods, Gaussian 
phase-space methods can be applied to either fermionic or 
bosonic quantum many -body svstems lfl6l fl7ll . These meth- 
ods employ a representation as a distribution over stochastic 
Green's functions. They have the essential property that they 
are probabilistic, allowing them to scale to large sized systems 
without the exponential complexity of the usual orthogonal 
basis-set methods. Such methods have proved useful for treat- 
ing the ground state of the fermionic Hubbard model lfl8l - l20ll . 
and dynamical cases where quantum Monte Carlo techniques 
are impractical J2ll l22tl . 

Yet one of the most fundamental properties of quantum sys- 
tems, entropy Jill, has not been treated using this approach. 
In ultra-cold atomic physics it is usually the entropy that is di- 
rectly measurable, rather than the temperature 112411 . This is 
because there is no traditional thermal reservoir, so that to 
measure thermal effects, it often necessary to use an entropy- 
conserving adiabatic passage to a state of known entropy. 
While some QMC methods can calculate entropy[|25], stan- 
dard phase-space simulations have not yet done this. Simi- 
larly, in quantum information, entropic concepts like entan- 
glement of formation [26] and the quantum discord lf27ll are 
important measures of quantum behavior. This leads to the 
question: can phase-space methods be used to calculate en- 
tropy, which is not a typical quantum observable? 

In this paper we investigate quantum entropy calculations 
using phase-space representations. Such calculations can be 
used, for example, to determine thermodynamics and entan- 
glement of a many -body system[28] or to calculate the fidelity 
of a quantum memory 12911 . Alternatively, they can be utilized 
to assess how close a calculated state is to a known state, or 
to check that unitary evolution preserves the entropy. This is 



essential in computationally demanding problems, where it is 
important to be able to check the validity of a given quantum 
simulation. More generally, such investigations throw light 
on one of the fundamental problems of quantum statistical 
physics. This is the well-known paradox that unitary evolu- 
tion leaves the quantum entropy unchanged, while apparently 
introducing disorder through collisions and mixing. 

For a probabilistically sampled distributions, we show that 
the preferred method for entropy calculations is a Gaussian 
phase-space representation. In this approach, the operator ba- 
sis is comprised of Gaussian operators, and the correspond- 
ing variance or Green's function is used as a stochastic phase- 
space variable lfl2tl . For ease of computation we treat the sim- 
plest entropy measure, the Renyi or linear entropy l30[l3ltl . An 
essential ingredient in sampled linear entropy calculations is a 
knowledge of the inner products of the generating operators of 
the phase space representation. We calculate these inner prod- 
ucts for phase-space representations of either bosonic [15]or 
fermionic lfl3ll fields. In the case of fermionic operators, the 
results come from an elegant application of Grassmann inte- 
gration in a Grassmann space. 

This method is directly useful for calculating a coarse- 
grained or localized entropy |[32tl . Coarse-grained entropy can 
increase even when the fundamental quantum entropy is in- 
variant, and is a fundamental entanglement measure. Coarse- 
graining is then carried out simply by restricting or project- 
ing the stochastic Greens function onto a subspace of sys- 
tem modes. As physical applications of these methods, we 
calculate the coarse-grained partial entropy for systems of ei- 
ther bosonic or fermionic modes that are coupled to a general, 
non-Markovian reservoir. These cases cannot be treated using 
master equation methods, and have applications to many cur- 
rent nanoscale systems in quantum information. These can be 
readily solved for the entanglement between the system and 
reservoir using the phase-space entropy approach. 

By comparison, the commonly used Wigner[i l]] and Husimi 
Q-function [2] expansions are not directly useful in entropy 
calculations. In the Q-function case, the inner products of 
the generating operators are not defined. In the Wigner dis- 
tribution case, the entropy obtained from a sampled phase- 
space distribution is singular. For normally-ordered distri- 
butions, like the Glauber-Sudarshan P-distribution(3l[33t], the 
sampled quantum entropy is well-behaved, but the distribu- 
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tion is singular, except for classical-like states. This rules 
out the use of traditional, classical phase-space mappings for 
entropy calculations using probabilistic sampling. The posi- 
tive P-distributionfl l34l 13511 . defined on a double-dimension 
phase-space, has a well-defined positive distribution with non- 
singular inner products. Even in this case, entropy calcula- 
tions are nontrivial since, for some quantum states, the sam- 
pling error for entropy calculations diverges. 

This paper is organized as follows: In the next section, we 
describe the definition of quantum entropy and also describe 
the phase-space representations that we will treat in this pa- 
per. We also give a general expression of entropy in terms of 
phase space-representations and discuss the evaluation of the 
sampled entropy. In section III we discuss coarse graining and 
reduced entropy and give their expressions in terms of phase- 
space projections. In section IV we describe sampled entropy 
calculations using fixed variance phase space methods, for in- 
stance the general Cahill-Glauber distribution, as well as the 
Husimi, Wigner, Glauber-Sudarshan and positive-P represen- 
tations. In section V and VI we discuss general Gaussian 
phase-space representations for bosons and fermions respec- 
tively. The evaluation of inner products of Gaussian operators 
is presented, which allows us to evaluate the linear entropy 
and the reduced coarse grained entropy. We evaluate the lin- 
ear entropy analytically for thermal states and a comparison 
with the sampled entropy using the Glauber-Sudarshan repre- 
sentation is made. Finally, as physically relevant examples, 
we evaluate the coarse grained entropy for both bosons and 
fermions, in the important case of a system of particles lin- 
early coupled to a non-Markovian reservoir. Section VII gives 
a summary of our results and conclusions. 



II. ENTROPY AND OPERATOR REPRESENTATIONS 

Entropy is a conserved quantity for unitar y ev olution in 
quantum mechanics. Intuitively, entropy 11231 13611 is a mea- 
sure of loss of information about a physical system. It is in- 
variant for unitary quantum evolution, and only changes when 
one considers a subsystem coupled to a reservoir. Physically, 
quantum pure state evolution involves no intrinsic information 
loss. Of course, this is somewhat counter-intuitive. One might 
expect the mixing effects of nonlinear evolution to reduce 
information. But this is only true if a restricted or 'coarse- 
grained' set of measurements is used; in principle, there is no 
information loss in pure state evolution, and information or 
entropy should be invariant in any simulation of unitary quan- 
tum dynamics. Thus, a pure state will remain a pure state, 
implying that state purity is a completely general benchmark 
for the accuracy of a quantum simulation. 



A. Quantum entropy 

Quantum entropy is most commonly defined using the von 
Neumann ll23l or Shannon entropy l37ll : 



This is an important physical quantity, related to both in- 
formation content and to thermodynamic behavior. However, 
there are many other conserved entropic quantities. The ex- 
istence of these can be thought of as related to general in- 
formation conservation. These are recently discussed [38] in 
relation to the Wigner function, where it is pointed out that 
any quantity like Sf = Tr (F(p)) is also conserved. In par- 
ticular, the linear entropy [39], which we normalize following 
Renyi[30], is defined as: 



S 2 = -ln7> (p 2 



(2.2) 



This has similar properties to the entropy, and measures state 
purity, since S2 = for a pure state, while S2 > for a mixed 
state. In this note we focus on the linear entropy, which is 
simplest to calculate using phase-space methods. This is also 
true for the fidelity of p to a fiduciary state po, 



(2.3) 



F = Tr(pp Q ) 



which is a closely related concept. Such fidelity measures [4C 
are useful in evaluating the accuracy of information storage in 
a quantum network l4lll . quantum computer !^ , or quantum 
memory B43II . 

The most general class of entropies normally studied in this 
way are the general Renyi entropies l30l [3lll . defined for p > 1 



S„ = 



1 



ln7>(p p ) . 



(2.4) 



It is known that S = limp_>i S p , so the conventional von Neu- 
mann or Shannon entropy can be regained from the general 
Renyi entropy in the appropriate limit. 



B. Phase-space representations 

A general phase-space representation can be written as[4, 



p = j P{X)A(X)dX 



(2.5) 



-TY(plnp) 



(2.1) 



where P(X) is the probability density over the phase-space, X 
is a real or complex vector parameter in a general phase-space, 
dX is the integration measure, and A(A) is the representation 
kernel or operator basis. For simplicity, we exclude phase- 
spaces that involve Grassmannian degrees offreedom ||45ll46Tl . 
These can be extremely useful in analytic calculations, but are 
not readily sampled computationally, since the vector param- 
eter X is not a real or complex vector. 

We consider a bosonic or fermionic quantum field the- 
ory with an M-dimensional set of mode operators a* = 

a\ , a\ , ■ ■ ■ .In the bosonic case, we can define 8a = a a 

and da' — a T — j3 ^ as operator displacements, where in gen- 
eral a and ^ are independent complex vectors. In the 
fermionic case we set these displacements to zero. The an- 
nihilation and creation operators satisfy (anti) commutation 
relations, with (+) for fermions and (— ) for bosons: 
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fl; a 



(2.6) 



The phase-space representations we will treat in this paper 
use a general number-conserving Gaussian operator basislll3l— 
[l5tl . in which any density matrix p is expanded in terms 
of a basis of Gaussian operators, defined as exponentials of 
quadratic operator forms A(A), where: 



A(A) = — A„ (A) 



1 



: exp 



8a f n8a 



(2.7) 



,yV"" y "' Jf 
Here, fl is a complex M xM matrix so that A = Ct,fi\fl 



,Jf = Tr 



A«(A) 



is a normalizing factor, and : : indicates 

normal ordering. The normalizing factor has two forms, for 
bosons and fermions respectively: 



det 
det 



-l 



2/ jU 



(2.8) 



The interpretation as a stochastic Green's function comes 
from the identification that fl is closely related to a correlation 

function of each basis member A(A): 



-T 



y 



21-11 



-T 



(2.9) 



In either case, the stochastic average of n over the distribu- 
tion P is physically a normally ordered many-body Green's 
function, so that: 



h & j) = { n ij + Pt a j)i 



(2.10) 



In traditional, classical types of phase-space - for example, 
the Wigner function approach - the random variable or phase- 
space coordinate is a stochastic position or momentum. In the 
case of a general Gaussian phase-space, the random variable 
is a stochastic correlation function. 

We see immediately, from Eq. ( 12.2b . that the Renyi entropy 
in a phase-space representation is: 

S 2 = -ln JJ P(X)P(X')Tr (a(X)A(X')^ dXdX' . (2.11) 

The evaluation of inner products of Gaussian operators of 
form Tr (a(X)A(X')\ is therefore a central task in calcula- 
tions of linear entropy using phase-space representations. 



C. Sampled entropy 

For computational purposes, distributions always exist such 
that P(X) has positive values, and it can be interpreted as a 



probability in these cases. One can then sample the distribu- 
tion N times to obtain a sampled estimate ps, such that: 



(2.12) 



This approximation becomes an exact equality in the limit of 
N — > oo, provided the sampling is unbiased. Given a set of 
samples A,-, we can now calculate the linear entropy as fol- 
lows: 



S2 



In 



(2.13) 



This, however, requires a double sampling of the popula- 
tion. In phase-space representations the kernel of the repre- 
sentation consists of non-orthogonal operators, so that the op- 
erator inner product Tr ^A(A;)A(A ; )^ is non-vanishing even 

if Xi ^ Xj. It is desirable that the two sets of samples A,, Ay 
are independent of each other, to prevent sampling biases. The 
above result has obvious extensions to other entropies. For ex- 
ample, the general Renyi entropy involves a p-fold summation 
over sampling indices: 



1 



■In 



^lTr(A(A yi )...A(A-; 
j=i 



,(2.14) 



but clearly the linear entropy is computationally the simplest. 
In the remainder of this paper, we focus on the question of how 
to evaluate the inner-products of the Gaussian phase-space ba- 
sis set, which is the essential ingredient in calculating a linear 
entropy or fidelity measure, and how to apply this in physi- 
cally relevant situations. We note that for some calculations it 
is useful to allow P(X) to have complex values B47I1 . in which 
cases the entropy is best calculated analytically. 



III. COARSE GRAINING AND REDUCED ENTROPY 

There is a fundamental paradox in understanding quantum 
entropy. For an isolated quantum system, all of the entropies 
defined above are invariant under unitary evolution, even in- 
cluding particle-particle interactions. This appears to defy 
conventional wisdom, which is that for a many-body system 
the effect of particle collisions is to cause mixing, and hence 
increase disorder. Thus, collisions appear to increase the en- 
tropy. Such expectations contradict the entropy invariance of 
unitary evolution, which is at the heart of such famous con- 
troversies as the Bekenstein-Hawking black hole information 
loss paradox SSI. 

However, these paradoxes are easily resolved at a practical 
level. Typically, in many experiments only part of the den- 
sity matrix is measurable. For example, one may only have 
experimental access to measurements of the low momentum 
modes. Under these conditions, one can separate the Hilbert 
space into a measured part Ma and unmeasured part Mb, so 
that the entire Hilbert space is M = M^ (£> Mb- Other sepa- 
rations of measured and unmeasured operations are also pos- 
sible, using the method of communication alphabets 115011 . An 
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interesting recent proposal of this type is to employ the many- 
body energy eigenstates as a communication alphabet to de- 
fine entropy 115111 . 

Here, for definiteness, we focus on a division of the Hilbert 
space into measured and unmeasured single-particle modes. 
These could, for example, correspond to a physical parti- 
tion into distinct spatial locations. The two parts of the 
quantum wave-function then become entangled during time- 
evolution under a Hamiltonian that couples the two parts. This 
means that part of the quantum information is only accessible 
through measurement of correlations. An estimate of this rel- 
ative entropy 15211 based on measurements reveals an apparent 
increase in entropy, or loss of information due to entanglement 

If we trace out the unmeasured part of Hilbert space, denot- 
ing this trace over j£g as 7>#, we obtain the reduced density 
matrix that corresponds to operational measurement on A: 



Pa = Tr B \ 



(3.1) 



Such a reduced density matrix can experience increased en- 
tropy - called entanglement entropy - even when the total en- 
tropy is conserved. The corresponding reduced entropy is 
then: 



vved 



1 



1 



; 1" ?>a(p£) 



(3.2) 



This reduced entropy is an important measure of quantum 
entanglement. In the case of a pure state, S'p d > is both 
necessary and sufficient for entanglement. This can also be 
extended to the case of mixed states. In this case, one must 
generalize the approach, to take account the possibility that 
the original state was a mixed state l26ll . 



A. Phase-space projections 

In the case of Gaussian phase-space expansions, all our en- 
tropy results are also applicable to the reduced entropy, in 
which case we must replace the phase-space basis A(A) by: 



A 4 (A ) = Trg A(A) 



(3.3) 



With such a replacement, the trace used in the following cal- 
culations must be replaced by a reduced trace over Ma for 
consistency. If coarse-graining is carried out on a modal 
basis, we can divide up the modes into two sets: a = 
[01,02, . ..%] = [a A ,a B ] . Here the modes cr may comprise 
only low-momentum modes, or alternatively, only modes lo- 
calized to part of an apparatus. 

In such cases, A A (X) depends on a new set of parameters 
k A = (ct A , fi A \ H A ) . The reduced displacements are just the 
projection of the full displacements onto the reduced Hilbert 
space, while the reduced covariance can be evaluated using 
standard trace identities. We first write the original matrix /i 
in a block form as: 



^AA pAB 



(3.4) 



so that the Gaussian exponent term becomes: 

8a^8a = da^fl^da* + So* ^ AB 8a B + {A^B) (3.5) 

Next, the relevant traces over the unobserved subspace B 
are evaluated using coherent state identities: 



1 



Tr b [0] = ^Jd™a(a\0\a), 
Tr f [d] = J d 2M a{-a\6\a) (3.6) 

This gives the result that the reduced basis set remains 
Gaussian, but with a modified covariance: 



A A (A)=A A (a A ,j3 At V) 



(3.7) 



where the reduced covariance matrix jtf A is given, for bosons 
and fermions respectively, by: 



BB 



BA 



21 jU 



BB 



BA 



(3.8) 



The important result here is that for a Gaussian basis, 
coarse-graining via mode-projection leaves the phase-space 
representation invariant. Just as for the full Gaussian expan- 
sion, there is a reduced Green's function for these Gaussian 
operators in the subspace. From Eq. ( |2.9t in the previous 
section, this must have the standard form of: 



(3.9) 



The results for the reduced stochastic Green's function for 
bosons and fermions can now be written, using standard ma- 
trix block reduction algebra, in terms of the block represen- 
tation of the original stochastic matrix n, which also has a 
decomposition: 



n A - 
Ltb — 




-T 


t 


n A - 


2/ A 


- E A 


-T 



n XA n AB 

n BA n BB 



(3.10) 



We find that the trace reduction simply gives the diagonal 
block in the first quadrant. 



A „AA 
l -f=H-f ■ 



(3.11) 



This has a simple physical explanation. We naturally expect 
that any correlation function that is restricted just to the A 
Hilbert space will have no dependence on measurable correla- 
tions of the B Hilbert space. This physical property of the full 
Green's function also holds for the stochastic Green's func- 
tions as well. 

The basis is mapped to new values, and the reduced density 
matrix p = 7>s[p] of Eq. ( 13.21 i can be written in the reduced 
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Gaussian representation in terms of the reduced Green's func- 
tion as: 



P A = Tr B J P{n)A{n)dn 
= j ' P(n)A A (n A )dn 



(3.12) 



Next, we can introduce the corresponding reduced distribu- 
tion function: 



P(n)dn 



(3.13) 



where n} A is the relative complement of n A , i.e., the set of all 
variables in n that are not included in « A . 

Using the definition of the reduced density matrix, the lin- 
ear coarse grained entropy, Eq. (13.21 ). for p = 2, is: 

Sf = -In ffp A (r^)P A {W A ) Tr A (a a (n A )A A (n ,A ) ) dn A dri A . 

(3.14) 

IV. FIXED VARIANCE PHASE-SPACES 

To evaluate the entropy from a set of phase-space samples, 
we need the inner-product of the kernel members. This de- 
pends on how the phase-space is parametrized, either through 
changing the displacement, or the variance, or both. Tradi- 
tional phase-spaces for bosons utilize a displacement-based 
approach, which is the most similar to classical phase-space 
ideas. In the case of fermions, the displacements must be 
Grassmann variables, not c-numbers 114511 . which means that 
only the variances can be readily sampled computationally. 
In this section, we treat fixed variance phase-spaces, which 
are therefore bosonic. 



A. Cahill-Glauber phase-space 

The traditional mappings of bosonic fields to a classical 
phase-space utilize a single classical displacement. These can 
all be written in a unified form as lfl5ll : 



A.v(A) = Jp:exp 



-(a t -a t )/i(o-a) 



(4.1) 



Here a,a^ are a complex vectors, and jU is held constant 
so that X = CC. There are three famous cases, correspond- 
ing to different values of fl = 2//(s + 1), where s = 0,±1. 
Cahill and Glauber 115 311 have calculated the inner product for 
these i-ordered representations, which includes the diagonal 
P-representation (s = 1) and the Wigner representation (s = 0), 
as special cases. Their results are that, for s > 0: 



7V(A,(a)A,(fl0 



1 

-exp 



a-a /s 



(4.2) 



We note that the Husimi representation with s = — 1 has no 
well-defined inner-product for its basis set members, as the 



product trace is divergent. Thus, a point-sampled Q-function 
is not a useful way to calculate the entropy, without additional 
assumptions. More sophisticated techniques would be needed 
in this case. One could, for example, expand the Q-function 
using Gaussian wavelets, instead of delta-functions, so that 
the sampling expansion uses smoother functions. However, 
since different types of s-ordering are interrelated through 
Gaussian convolutions, this simply generates another member 
of the class of Gaussian operator expansions. 



B. Wigner representation 

Strictly speaking, the only positive Wigner distributions are 
the Gaussian ones that represent certain special cases, includ- 
ing thermal, coherent and squeezed states. Nevertheless, one 
often wishes to use a truncated Wigner time-evolution equa- 
tion, which generates positive Wigner distributions as an ap- 
proximation to the full time-evolution. This has a close anal- 
ogy with a classical phase-space, for which entropy can also 
be calculated in the classical sense. 

One can treat the Wigner case as the Cahill-Glauber repre- 
sentation in the limit of s — > 0, where: 



Tr Ao(a)Ao(a') =n M 8 M (a-a 



(4.3) 



which is highly singular. As in the Husimi case, point- 
sampling doesn't provide a useful estimate of the purity. Two 
distinct samples will not have identical points in phase-space, 
except for points of measure zero where the samples are equal. 
This demonstrates the nontrivial nature of estimating quantum 
entropy in sampled phase-space representations. One can un- 
derstand this from the perspective that the coherent states are 
the only pure states with a positive Wigner function. These 
have a finite distribution variance, but a zero quantum entropy. 

This result is consistent with other calculations. It is known 
that one can estimate in a Wigner representation through 
S: 

S 2 = -\n7t M Jw 2 (a)d 2M a 1 (4.4) 

which is identical with the delta-correlated trace result given 
above. This can be used when W is known analytically, but it 
is not computationally useful when we only have access to a 
sampled estimate of W (a). 



C. Glauber-Sudarshan 

In the case of the normally-ordered Glauber-Sudarshan rep- 
resentation, s = 1. This corresponds to an expansion in 
coherent-state projectors, so that Ai(flt) = |0t) (<x\, where |ot) 
is a coherent state, and 

(4.5) 

Here the linear entropy is well-behaved, and both linear 
coupling and damping can be treated exactly. However, there 
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Figure 1. Comparison of the linear entropy for thermal states using 
the Glauber-Sudarshan representation and the Gaussian representa- 
tion for bosons. The dotted line is the exact result using the Gaussian 
representation for bosons for n = 0.1, 1, 10, 100, 1000 , with smallest 
occupation numbers having the lowest entropy. Here N is the number 
of samples used. 



is no corresponding stochastic process in this case, for non- 
linear evolution of an interacting system, and many nonclas- 
sical states involve a nonpositive or singular distribution. For 
a positive distribution, the only pure states in this representa- 
tion are coherent states. Provided a Glauber-Sudarshan dis- 
tribution exists, a direct point-sampling is enough to obtain 
the entropy. One can easily obtain the entropy of a non- 
interacting thermal state, which always has a well-defined 
Glauber-Sudarshan distribution. For example, the vacuum 
state has a delta-function distribution, and so clearly one has: 

Tr ^Ai (tt)Ai {ct'^j = 1 , and hence S2 =S P =0 as expected. 

For a thermal case with: 



Pth x exp 



exp 



-n] 1 a 



(4.6) 



where n is the thermal Bose-Einstein occupation number, 
clearly 



e E k /k b T _ 1 



, where kj, is the Boltzmann constant. 

Here one finds that in the Glauber-Sudarshan representa- 
tion, one has: 



P(a) = exp [• 



a f « 'a' 



(4.8) 



Therefore we can use the results of P(ot) in order to sam- 
ple the entropy for the thermal states. In Figure ([TJ we show 
the results of the sampled linear entropy as a function of the 
number of samples N for different values of the thermal Bose- 
Einstein occupation number n, compared with exact results 
obtained in the next section. 

The generators in this case are coherent state projection op- 
erators, which means that obtaining a coarse-grained entropy 
is straightforward. On dividing the modes into two groups, 
A and B, one can simply write the coherent state as an outer 



product: 



so that: 



\a) = \a A )\a B ) 



A 1 (a)=A 1 (a A )A 1 (a B ) 

where Ai(a A ) = \ct A ) (& A \ is a coherent state projector in the 
reduced Hilbert space. 



D. Positive-P representation 

The positive P-representation extends the Glauber- 
Sudarshan representation into a space of double the classical 
dimension, with A — ( a , P ) . This the advantage that any state 
or density matrix has a positive probability expansion. Unlike 
the Husimi Q-function, the basis set has a non-singular 
inner-product, which allows the entropy to be calculated 
through sampling techniques. The kernel can be written in an 
alternate form as an hermitian projection operator 0] : 



A(A) = 



\*)<fi\ 

(P \ot) ■ 



(4.9) 



Just as in the Glauber-Sudarshan case, the issue of coarse- 
graining is a straightforward one of simply dividing the modes 
into two groups, and replacing A(A) by its reduced version, 
A(A ). The inner-product is always well-defined, being just a 
Gaussian form in the displacement vectors: 



Tr ( A(A)A(A'' 



(fi\a')(fi'\a) 
(P\a)(p'\a>) 



exp 



(a a 



.(4.10) 



In all cases, a highly localized distribution is guaranteed to ex- 
ist, from the fundamental existence theorem of the positive-P 
representation. This states that at least one canonical, positive 
distribution P(a,j3) always exists for any p, with: 



(4.7) P(<*,P) 



1 



(2n) 



2M' 



-\a-P\ z /4 



a + p 



a + p 



(4.11) 



While this distribution always exists, and is suitable for cal- 
culating moments, it generally leads to large sampling errors 
when calculating the entropy. This is due to the fact that when 
I Of; — Pi\ 2 3> 1, in Eq. ( 14.10b . the cross-terms can become ex- 
ponentially large, as these are not sufficiently bounded by the 
exponentials in the canonical form, Eq. ( 14.1 It - 

In summary, we see that calculating entropies using a dis- 
placement based phase-space expansion with point sampling 
is non-trivial. With traditional phase-space expansions, ei- 
ther the basis has singular inner-products, or the distribution 
is non-positive, or both. For the Glauber-Sudarshan represen- 
tation of a thermal state, the distribution is well-behaved and 
the linear entropy can be computed. In the positive-P case, a 
positive distribution always exists, and the basis has nonsin- 
gular inner products. However, even in this case the entropic 
sampling error can diverge for nonclassical states. 
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V. GAUSSIAN REPRESENTATIONS FOR BOSONS 

An alternative way to represent quantum states in phase- 
space, is to use a general representation in terms of Gaussian 
operators. These types of phase-spaces can in principle com- 
bine the displacement and variance-based approach. How- 
ever, for definiteness, in this section we will treat the case 
where the representation is entirely variance based. Such an 
approach has a clear intuitive meaning. In this approach, the 
physical many-body system is treated as a distribution over 
stochastic Green's functions, whose average in the observed 
Green's function or correlation function. We note that the ba- 
sis set includes non-hermitian terms for completeness, which 
means that the stochastic Green's functions themselves can be 
nonhermitian. 



A. Un-normalized Gaussians 

For the bosonic case, we must evaluate the trace of the 
product of two un-normalized bosonic Gaussian operators, 



B 



(W 



Tr 



A« M A » (V) 



for the M-mode case: 



B[H,V) = Tr 



,-a fia 



(5.1) 



Using the expressions for the trace of an operator, Eq. ( 13.61 ). 
and the expansion of the identity operator in terms of the 
bosonic coherent states: 



±l<l^a\a <a\=l. 



(5.2) 



we obtain: 



x\p)(P\:e-^:\a). (5.3) 



Expanding the normal-ordered exponential, and using the 
standard eigenvalue properties for the bosonic coherent states: 
a a = a\a) gives: 



B = ^fd™ ad™ ft {a\e-"' 
|/3)</3| e -^«|a>. 



(5.4) 



From the inner product of two coherent states, we finally ob- 
tain a Gaussian integral over 2M complex coordinates: 



t ^3-/J t va-|a-/J|- (55) 

Next, introducing a double-dimension vector: 

a 

7 = P 



(5.6) 



we can write this as: 



det[r; 



(5.7) 



where we have used the standard identity(||54j)) for an N- 
dimensional Gaussian complex integrals, and introduced a 
double-dimension matrix, 



l m-z: 

v / / 



(5.8) 



Therefore, on simplifying the determinant, we obtain: 

-l 



B /i,V = det 



I- M-/ (v-L) 



(5.9) 



B. Normalized Gaussians 

It is useful to rewrite these expressions in terms of the un- 
derlying stochastic Green's functions. These are the normally 
ordered correlations of the basis sets, defined so that: 



nu = Tr 



(5.10) 



Using this definition, the normalized Gaussian generators are: 
1 



A(n) 



det[I + n] 



: exp 



(5.11) 



We note that there is a restriction on the values of n , which 
is that 3i {I + n} must have positive definite eigenvalues in or- 
der for the basis operators to be normalizable, and hence for 
the Gaussian generators to be in the Hilbert space. 

Applying this normalization to the results given above, one 
finds that: 



Tr A(n)A(m) = det [I + n + m] 



-l 



(5.12) 



For Renyi entropy calculations there is another restriction. 
This is that all pairs of stochastic samples must have the prop- 
erty that 9J{I + n + m} has positive definite eigenvalues to 
calculate the entropy using sampling methods. Under this re- 
striction, the inner products are well-defined. 

In order to illustrate the technique of the Gaussian repre- 
sentation for bosons in the evaluation of the linear entropy 
and coarse grained entropy respectively, in the next two sub- 
sections we will evaluate the linear entropy of thermal states 
and the coarse grained entropy of a system coupled to a non- 
Markovian reservoir. 



C. Thermal linear entropy for bosons 

The linear entropy for thermal states is evaluated as pre- 
viously, using Eq. ( 12.13b and the result of Eq. ( 15.121 ) for the 
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Si 


0.01 


0.0198 


0.1 


0.1823 


1 


1.0986 


10 


3.0445 


100 


5.3033 


1000 


7.6014 



Table I. Linear entropy for thermal states using the Gaussian repre- 
sentation for bosons. 



single-mode case. When the density matrix is thermal, only a 
single basis set member is required, and: 



S 2 = -\nTr (A> rft ) 
= rndetg + 2n, ft ] 



(5.13) 



For the single-mode case we know that the thermal Green's 
function is a scalar: ry h = n t h, where is the Bose-Einstein 
occupation number at a given temperature. In Table 1 we show 
the results for the linear entropy using the Gaussian represen- 
tation for bosons. In Figure (fl} we show the comparison of 
the results using the Gaussian phase-space representation, re- 
sults of Table(Q]i, and the results from the sampling using the 
Glauber-Sudarshan representation as a function of the number 
of samples, giving excellent agreement in the limit of large 
numbers of samples. 

Clearly there is a great improvement in efficiency in this 
case, relative to the Glauber-Sudarshan approach. Only one 
Gaussian phase-space sample is needed, instead of up to 
50000 samples using more traditional phase-space methods. 



where 03% = d&Vk + gik- Here vjt denotes the resonant fre- 
quencies of the modes, g^i denotes the couplings between the 
modes. We assume that the modes for k = l,...S are system 
modes (A), while the remainder are the reservoir (B). We note 
that we make no assumptions concerning their relative sizes 
or quantum states. We suppose that the initial density matrix 
at time t = has the general number-conserving form: 



P°(n)A(n)dn 



(5.15) 



The Gaussian representation provides a form to express the 
real or imaginary time evolution of the density matrix of ei- 
ther fermions and bosons into a set of phase-space stochastic 
equations HUH EH, the mappings are given by lfl4l[T5ll : 



n ik p 
pn ik 



"at 



dn b 
d 



diti, 



-(1 ± «;,„)«« 
-n im (l ±n a ) 



(5.16) 



Here the +(-) sign corresponds to the bosonic (fermionic) case 
respectively, and we will use the bosonic identities here. 

The real time evolution of the density matrix is given, as 
usual, by: 



dp 
~dt 



(5.17) 



Using Eq. ( 15.171 ) and Eq. (15.16b we obtain that the time evolu- 
tion equation of the stochastic Green's function in matrix form 
is: 



D. Coarse-grained entropy for bosons 

We now wish to consider a practical example of consid- 
erable physical applicability. In much of modern physics a 
bosonic mode is coupled to a reservoir, with which it can ex- 
change particles. The traditional example of this is a single- 
mode interferometer l55ll . In current applications relevant to 
quantum information, one may have a localized photonic 
waveguide mode 1 56], an ultra-cold Bose condensate lf57ll . or 
a nano-mechanical oscillator ll58ll . These exchange photons, 
atoms or phonons respectively with their environments. In 
nearly all of these recent applications, one is interested in evo- 
lution with non-Markovian reservoirs. 

In order to model such physically important examples and 
evaluate the coarse grained entropy, we consider the following 
non-Markovian system. A set of bosonic modes (the system) 
is coupled to a large number of other modes (the reservoir). 
The total Hamiltonian can be written, on introducing hn = 
ajdj, as: 



» = (5-18) 
The solution of Eq. ( 15.181 ) is simply: 

n(t) = e im n(0)e- im . (5.19) 

The reduced linear entropy for this system is evaluated using 
Eq. (13.141 ) and the result for the time evolution of the stochas- 
tic Green's functions, Eq. ( 15.191 ). In this case we will trace 
over the system A. Therefore, the expression for the reduced 
entropy is: 

S r f = -In jj P°(«)P°(n')det[l +n A (t) +^ A (t)]- l dndti. 

(5.20) 

E. Example of bosonic entropy 



ik 



(5.14) 



In order to illustrate the time evolution of the coarse grained 
entropy of Eq. ( 15.201 ). we will consider the following model 
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1.4r 

1.2 




50 



couplings needed to obtain an exactly soluble form, there are 
no restrictions to the state, the type of coupling or the subdivi- 
sion between the system and the reservoir in this calculation. 

VI. GAUSSIAN REPRESENTATIONS FOR FERMIONS 

The fermionic case is similar, except that one must use 
fermionic coherent states [45] and Grassmann integrals to 
carry out the trace calculations. Just as with bosons, this has a 
clear intuitive meaning. In this approach, the physical many- 
body system is treated as a distribution over fermionic Green's 
functions, whose average in the observed Green's function or 
correlation function. As with the bosonic case, the stochastic 
Green's functions themselves can be nonhermitian. 



Figure 2. Time evolution of the coarse grained entropy for a bosonic 
thermal state coupled to a non-Markovian reservoir. 



for the Hamiltonian of Eq. d5.14l >. The bosonic system will be 
a single mode (j — 1) and the reservoir will be modeled as a 
Lorentzian distribution of couplings, with: 



8i 



(5.21) 



where C is the strength of the coupling, v, = ±jdco are the 
resonant frequencies of the modes, and s describes the non- 
Markovian reservoir width. For this model, the time evolution 
of the stochastic Green's function of Eq. ( 15.191 ) is written as- 
suming that «(0) describes the thermal state at t = 0, with a 
system occupation of n = 1 , and all the other modes unoccu- 
pied. Here ffl is an M x M-matrix, and M is total the number 
of modes, so that: 



/ Vn g 2 \ glM 

812 '■■ ••• : 



0) 











•. 



\ 8MI 



(5.22) 



VMM ) 



In Figure (0, we show the time evolution of the coarse 
grained entropy. We use the following parameters: M = 100, 
s = 0.5, C = 0.05. The frequency spacing of the modes is 
dco = 0.04. We observe that the non-Markovian behavior of 
the reservoir, as shown in the increase and decrease of the en- 
tropy with time. 

In summary, to solve for the coarse grained entropy one 
must take the partial determinant of a block reduced form n A 
of the time-evolved stochastic Green's function, then average 
over the initial ensemble. Apart from the limitation to linear 



A. Un-normalized Gaussians 

Here we must evaluate the trace of the product of two 
un-normalized fermionic Gaussian operators, ^^jU,V^ = 

Tr A u (ft) k u (v) for the M-mode case: 



(M,v) 



Tr 



,-a< fia .. -a va 



(6.1) 



For fermions ll45Tl . the trace of an operator using fermionic 
coherent states la) in terms of Grassmann variables a is: 



Tr[0] = j d 2M a{-a\6\a) 1 

and the identity operator is: 



(6.2) 



d ZM a\a)(a\ = l. 

Therefore, we have: 

x\P){fi\:e-*#:\a). 



(6.3) 



(6.4) 



Expanding the normal-ordered exponential, and using the 
standard eigenvalue properties for the fermionic coherent 
states: a\a) = Ct\a) gives: 



f (m,y) - J d 2M ad 2M p (-a\e at aP 



\P)(P\ 



,-/3 f vo 



a). 



(6.5) 



From the inner product of two fermion coherent states, we 
note that: 



(a\P)=e atp -( attt+p ' p ) /2 



(6.6) 
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Next, introducing a double-dimension Grassmann vector : 



(6.7) 



we finally obtain a Gaussian Grassmann integral over 2M 
complex coordinates, which we can write this as: 



j4My e « 1 M^-^ 1 va-a t /3+^ 1 a 1 -(a 1 a+J3 t J3) 



= det[r] (6.8) 

Here we have used the standard identity([54]) for an N- 
dimensional Gaussian complex Grassmann integrals, and in- 
troduced a double-dimension matrix, 



L Z-/£ 
V— 7 I 

Therefore, on simplifying the determinant, we obtain: 



(6.9) 



F fl.V 



det 



i+ 7-M )(L-v) 



(6.10) 



B. Normalized Gaussians 

Just as in the bosonic case, it is useful to rewrite these ex- 
pressions in terms of the normally ordered Green's functions 
or correlations of the basis sets, defined so that: 



mj = Tr 



A(n)a ( 'fl ; - 



Here, introducing the hole Green's functions, ft = [I — n], and 
10 = [I — m], the normalized generators are: 



A(n) = i2det [ft] : exp a [n -1 - 21] a 
and therefore, 



Hence, 



v T = 2i-fr 1 
H T =2I-m" 1 



fi(7-V) =n 



(6.12) 



(6.13) 
(6.14) 



(6.15) 



which leads to the following result for the normalized inner 
product: 



Tr 



A(m)A(n) = det [nm + nm] 



(6.16) 



We note that this has some obvious properties. Suppose 
that n and m are each diagonal in the same basis, with real 
eigenvalues n, such that < n; < 1 . Then one obtains: 



Tr 



M 



A(m)A(n) = ]~[(n,-m / + n,-m i ) (6.17) 



"th 


s 2 


0.01 


0.02 


0.1 


0.1984 


0.5 


0.6931 


0.9 


0.1984 


0.99 


0.02 



Table II. Linear entropy for thermal states using the Gaussian repre- 
sentation for fermions. 



Thus, the two generators are orthogonal if, in any mode, one 
generator has a vanishing particle population while the other 
has a vanishing hole population. The overlap is maximized 
if the generators both have a unit hole population or a unit 
particle population in all modes. 

For the thermal case, the entropy can be evaluated in other 
ways, but here we demonstrate the technique using the Gaus- 
sian operator method, which will be useful to evaluate the en- 
tropy of other systems. 



C. Thermal linear entropy for fermions 

We can now apply these inner products to the evaluation of 
the linear entropy of a thermal Fermi-Dirac states. When the 
density matrix is thermal, only a single basis set member is 
required, and: 



S 2 = -\nTr(A 2 (n th )' j 
= -lndet[l-2n, ft +2n, 2 J 



(6.18) 

(6.11) j us t as with bosons, for the thermal case we know that the 



thermal Green's function is a scalar: = «,/,, where n,/, is 
now the Fermi-Dirac occupation number at a given temper- 
ature, so that < rith < 1. Here the results are asymptoti- 
cally equal to the bosonic case as expected for n,/, C 1 or for 
nth -C 1 • Typical results are shown in Table (II), showing the 
particle-hole symmetry. The greatest entropy is at n t f, = 0.5, 
corresponding to infinite temperature, while mirror states with 
small hole occupations can be thought of as having negative 
temperatures or negative Hamiltonians. 



D. Coarse grained entropy for fermions. 

Just as in the case of bosons, we now consider an exam- 
ple of physical applicability for the case of fermions that is a 
fermionic mode coupled to a reservoir. An example of such a 
system is a quantum dot coupled to a fermionic reservoir 15911 . 
or a fermionic atom-tronic circuit! 60]. This system can be 
considered as an example of solid-state quantum physics and 
has potential applications in quantum information processing. 
In such hybrid quantum systems, long-range interactions can 
be important. Here we neglect this in order to obtain analytic 
results, although these can be added if necessary. 

We consider an identical Hamiltonian to the bosonic case 
in the last section. Similar to the bosonic case, we assume 
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that the modes for k = l,...S are system modes (A), while 
the remainder are for the reservoir (B) and the initial density 
matrix at time t — has the general number-conserving of 
Eq. EE). 

Using the identities of Eq. ( 15.161 ) and Eq. ( 15.17b we obtain 
the time evolution equation of the stochastic Green's function 
is identical to the bosonic case. Therefore, the expression for 
the reduced linear entropy for fermions is: 

S r f - - In jj P° ( n ) I* ' ( n' ) det \n' ( t )A ,A (t)+r£(t )n A (t )] Jndn! . 

(6.19) 

This has a very simple physical interpretation. The linear 
entropy and hence the fermionic entanglement can be calcu- 
lated completely from the local stochastic Green's functions 
in the system of interest. However, one must average over all 
possible initial states defined by the complete initial phase- 
space distribution P°{n). This is necessary, since the corre- 
lations and initial states of the reservoir can change the final 
system properties. 

E. Example of fermionic entanglement 

Similar to the bosonic case, we will illustrate the time evo- 
lution of the coarse grained entropy of Eq. ( 16.191 ). We will 
consider a pure fermionic number state. The (B matrix of time 
evolution of the stochastic Green's function has the same form 
as the bosonic case, Eq. ( 15.221 ). We also model the reservoir 
with a Lorentzian distribution described by Eq. ( 15.211 ). 

In Figure (0), we show the time evolution of the coarse 
grained entropy, Eq. ( 16.191 ), using an identical model and pa- 
rameters to the bosonic case. However, there is a large physi- 
cal difference, as a fermionic state with n = 1 is a pure num- 
ber state. We observe an initial increase of entropy, which 
means that the entropy is measuring the entanglement of the 
system with the reservoir. As before, the non-Markovian be- 
havior of the system is shown in the increase and decrease of 
the coarse-grained entropy, which in this case corresponds to 
entanglement oscillations. 

VII. SUMMARY 

We have calculated the linear entropy for sampled phase- 
space representations of bosonic and fermionic quantum 
many-body systems. The crucial element to the calculation 
is an evaluation of the inner products of the phase-space basis 
elements. Traditional displacement-based phase-space meth- 
ods have a range of pathologies. In the Wigner and Husimi 
cases, the inner products are singular or divergent, while in 



the Glauber-Sudarshan case, the representation is not well- 
defined in all cases. Even the positive-P distribution, which 
exists and has well-defined inner products, we find there is 
a sampling convergence problem. By comparison, Gaus- 
sian phase-space representations for fermions and bosons are 
much more suitable for the task. For thermal states, only a sin- 




10 20 30 40 50 

t 



Figure 3. Time evolution of the coarse grained entropy using the 
fermionic Gaussian representation, for a pure fermionic number state 
coupled to a non-Markovian reservoir. 



gle basis element is needed, and the inner-products are well- 
behaved. 

There is a counter-intuitive element to the idea that entropy 
is conserved in quantum dynamics; but this must be the case 
when simulating time-reversible, unitary quantum dynamics. 
We show how, in the case of reduced entropy of a subsystem, 
the linear entropy can and does evolve in time. We give an 
exact calculation of couplings of Fermi and Bose systems to 
a non-Markovian quantum reservoir. Such phase-space meth- 
ods appear useful for investigating the fundamental paradox 
of entropy-invariance in unitary quantum dynamics. They are 
equally applicable to entire system evolution and to the evolu- 
tion of the density matrix for a coarse-grained sub-space. 
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